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In the present talk I discuss the recent progress in the theoretical descrip- 
tion of the Deep Inelastic Scattering (DIS) at small x. From the theoretical 
point of view, the inclusive cross-section of DIS can be regarded as a prod- 
uct of the leptonic and the hadronic tensors, with exchange by a deeply 
off-shell virtual photon. Further on I discuss only the hadronic tensor. The 
hadronic tensor consists of two parts: the spin- dependent part and the spin- 
independent one. My talk is about the spin-independent part, W^^^. This 
tensor should be made of on-shell momenta P (the initial/final hadron mo- 
mentum) and of the deeply off-shell momentum q (the incoming/outgoing 
virtual photon momentum). Also it should bear the Lorentz indices n and 
u corresponding to the virtual photon polarizations. It must also respect 
both the gauge and the Lorentz invariance. This leaves us with two inde- 
pendent projection operators, each multiplied by a scalar function. These 
are the spin-independent structure functions Fj^ and F2. Their arguments 
must be some scalars made of P and q. Traditionally, they are chosen as 
g2 = _g2 > and X = gV2Pg, (0 < X < 1). Thus, 

W,. = (-g^^+^)F,\X,Q') + {P,-q,^){P^-q^^)^F,\X,Q') . (1) 
qz qz qz pq 

The next step is exploiting the factorization. According to it, the whole 
hadronic tensor, and its spin-independent part in particular, can be regarded 
as a convolution of two objects: the probability to find a parton (a quark 
and or a gluon) and the Deep Inelastic Scattering off the parton. Such par- 
tons are supposed to be nearly on-shell, so one can construct the "partonic" 
tensor for them in the same way as Eq. (0) was made, replacing P by p and 
Xhy X. Such new spin-independent tensor w^y (parameterized by "partonic" 
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structure functions Fi and F2, each of them depends on Q"^ and x = Q'^/2pq), 
can be calculated with the perturbative QCD methods and must be convo- 
luted with . As for calculating w^^, it is a purely technical problem, 
though quite complicated. However, it is scarcely possible to convolute w^^, 
with because presently there is no model-independent way to calculate 
Pj^. Instead, the regular way to obtain Fj^^ is imitating this procedure by 
picking up good fits to fill in the gap between W^^^ and w^y. From now on I 
discuss w^i,, i.e. DIS off a parton. I choose a quark as this parton. In order 
to calculate Wf^,y, one have to calculate all involved Feynman graphs and sum 
up their contributions. In the Born approximation 

Fi = (eJ/2)5(l-x), F2 = 2xFi. (2) 
where Cq is the electric charge of the initial quark. 

The radiative corrections to the Born approximation can be calculated in the 
flavour- dependent way where the quark interacting with the virtual photon 
is the initial quark and in the flavour-independent way where these quarks 
are different. In this case, there is no dependence on the flavour of the initial 
quark because one can sum over flavours of the quark interacting with the 
virtual photon. Therefore, any structure function has both flavour- dependent 
(non-singlet) and flavour-independent (singlet) parts. Basically, the singlet 
parts of the spin-independent structure functions are more important because 
they define the total cross-sections. However, the non-singlet parts are also 
interesting to know. For example, one can subtract a structure function for 
DIS off proton from the same structure function for DIS off neutron. The 
result is flavor- dependent. Then, the non-singlet structure functions and the 
singlet ones have some common features. On the other hand, the non-singlet 
structure functions are simpler objects to calculate. So, they can be regarded 
as a test field for checking various assumptions or models made for the singlet 
structure functions. The most important radiative corrections to the Born 
values (0) are logarithms of x and Q^. When x is large, logarithms of x 
are not important but when x is small, none of these kinds of logarithmic 
contributions must be neglected a priori. The less is value of x, the more im- 
portant are contributions of higher-order radiative corrections. With certain 
approximations, there were developed regular methods to calculate them to 
all orders in a^. The most wide-spread methods among them are the DGLAP 
and the BFKL 0. Both of them are one-dimensional linear differential 
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equations with integral operators in their right-hand sides. However, they 
were developed for operating in different kinematical regions: the DGLAP 
was suggested for studying the kinematical region a; ~ 1 and large . It 
neglects therefore logarithms of x which are not multiplied by logarithms of 
. So, it respects Q^-evolution only. For evolution of the singlet part of 
F2 = F the DGLAP equation can be written as 

dF 

PF (3) 



ainQ2 

where P is the DGLAP integral operator, with the splitting functions as a 
kernel. The splitting functions are known to the a1 -order. After the Mellin 
transform 

/ICO flfj , , , 

^JLe^^<V^)F{N^Q% (4) 
-100 

the DGLAP is 

dF{N, t)/dt = -f{N)F{N, t) (5) 

where t = InQ^. 7(A^) are the anomalous dimensions. They are known to 
the first two orders in as- 

The BFKL was done for accounting for x-evolution. It sums up Inx to all 
orders in a^- Like the DGLAP, it is also one- dimensional evolution equation: 

where K is the BFKL integral operator, and Fq is the singlet F2 in the lowest 
order. K is presently known to the order a^. After the Mellin transform (^), 
the BFKL for the gluon-gluon elastic forward scattering is 

NF = 5^^\h - k2) + {Ko + K^)F, (7) 

where Kq is the leading order (LO) integral operator and Ki is the non- 
leading order (NLO) one. In particular, is 



iin-T - 



TT J {ki - kf 



(8) 



for forward scattering. ki^2 in Eq- (0) are transversal (two-dimensional) ex- 
ternal momenta of the upper and the lower ladder gluons respectively. 
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Along with their advantages, both the DGLAP and BFKL have certain dis- 
advantages. The DGLAP is practically free of disadvantages when it operates 
in the region of large x but when it is extrapolated into the region of small 
X, contributions ~ In x'^ systematically neglected in the DGLAP to all orders 
in as must be incorporated. Without them, the standard DGLAP prediction 
for the asymptotically small-x behavior is 



exp 



4n, 



In 1 /x In 



ln(/iVA|cD) 



1/2 



(9) 



where ric = 3, b = (llric — 2ra/)/127r, Uf is the number of involved quark 
flavours and /i^ is a starting point for -evolution. 

The BFKL accounts for contributions ~ In*^ x to all powers in the QCD 
coupling. However from the theoretical point of view, the main disadvantage 
of the BFKL is that it treats as a constant. Thus, its small-x prediction 



F r^X-^^^ 



with the Pomeron intercept 

Ap 



TT 



41n2 



1 + r- 



TT 



(10) 



(11) 



cannot be used without specifying a value of as- In the MS'-scheme used in 
Eq. (pTD, r = —20 for four involved quark flavors and = 0.2. Thus, the 
NLO-corrections (the second term in (|Tl|)) are enormously big. Basically, 
recent works on the unpolarized DIS structure functions are intendant to 
overcome these disadvantages of the DGLAP and the BFKL with different 
means. 

The first possible solution to the problem of the big NLO corrections in 
Eq. (pAl) is suggested in 0. In difference to the MS- scheme used for ob- 
taining (|l3), they suggest to use the approach of work 0] on the optimal 
scale setting in (pJ]). With it, the value of as, still considered as fixed, is 
decreased down to as = 0.15 and at the same time the value of r becomes 
lesser: f ^ r/2. Another improvement of the BFKL in the framework of this 
approach is that the new scale does not lead to negative cross-sections for 
non-leading modes as it was when MS was used. 

On the contrary, as is supposed to be really running in works 0. The main 
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assumption made there is that as in the BFKL -kernel should depend 
on virtuality k\ of the external ladder gluon. The main conclusion is that 
when running as{k\) is taken into account, this new version of the BFKL 
predicts a small-x behavior of the DGLAP type rather then the power-like 
(the Regge-like) behaviour (p!0|). 

The QCD coupling in the BFKL- kernel is supposed to be running also in 
works 0, though its argument is more complicated compared to Indeed, 
Eq. (ID reads that the integration region over k includes the subregion (i) 
where 

k^ < kl (12) 
as well as the opposite subregion (ii) where 

e > kl (13) 

At last, there is a subregion (iii) where 

k^ ^ kl (14) 

Applying the Mellin transform (^ to the BFKL and, after that, applying to 
the result the new Mellin transform 

/lOO rl ]\/T 
—e^'F{M) (15) 

converts the integral equation (|^) into an algebraic one: 

NF = Fo + x{M)F. (16) 

x(M) in (0) is the characteristic function of the BFKL. As the integral 
operator K in (P) consists of the leading order contribution Kq and the non- 
leading one Kl, X = Xo + Xi respectively. The LO characteristic function 

n Nr 

xo = ^ mi) - m) - - M)] , (17) 

TT 



with ijj being the Euler F -function derivative. Eq. (|T^ reads that xo is 
invariant to replacing M by 1 — M, also it has singularities at M = and 
M = 1. With good accuracy. 



Xo ~ Xo 



col OisNc 



IT 



1 1 

M ^ 1 -M, 



(18) 
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Corespondingly, Kq can be approximated by K^"'' through the inverse Melhn 



transform of (|T^). This approximation leaves out kinematical region (14) 
but it is supposed in ^ to be not important. Suggesting that in remaining 
regions (0) and ([T3|) in Kq"'- should depend on the largest gluon virtuality, 
works assume that depends on kf in the region (0) and on fc^ in the 
region (|13[). Then, using similar arguments, works approximate the NLO 
BFKL characteristic function xi by 



Xf = (a.iVe/7r)^(-ll/12) 



1 1 

+ 



M2 (1-M)2 



(19) 



so that discrepancy between xi and XT'" ^^^^ than 7 per cents. The number 
(-11/12) in (|l^) is the subleading contribution to the leading order DGLAP 
gluon anomalous dimension = l/uj — 11/12. Approximations (|18|), (|1^) 
in works and their treat of ag are expected to imitate next-to- next order 
contributions (and so on) to the BFKL in a rather technically simple way 
and still to be in a reasonable accordance with the BFKL. 
The approach using the BFKL for improving knowledge of the DGLAP 
anomalous dimensions and splitting functions is suggested in works . The 
point is that, instead of straightforward procedure of accounting for non 
leading order contributions to jgg to all orders in as, one can use for it the 
BFKL which already includes such a resummation. After the double Mellin 
transform both with respect to x and with respect to In (we remind that 
we keep notations ^ = ln(l/a;) and t = lng^//i^ through this paper), 

/loo rl ]\/T ri'<^ rl M 
|^F(M)e*^*+^^F(M,Ar) (20) 
-too Ztti J -too Zm 

the BFKL is (cf. Eq. (|16D) 

NF{M, N) = Fo + x(M)F(M, A^) (21) 
with the obvious solution 

F(M, N) = — -. (22) 

Eq. ( p2D reads that there is the relation between singularities of F in M and 
in A^. Indeed, 

N = x{M). (23) 
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at points of singularity of Eq. (|2^). 

On the other hand, the DGLAP reads the asymptotically small-x behavior as 
F ~ exp[7(iV)t] at such points, which being compared to (^) immediately 
gives 

M = -i{N). (24) 
Combining Eqs. (|23|) and (p4|) leads to 



iV = x(7(iV)), M = 7(x(M)) (25) 



Obtained in [H, |^ Eqs. ( pS] ) are called the duality relations. As x contains 
contributions ~ {ag/M)^ to all powers in the relation iV = xili^)) can be 
used to express the anomalous dimension 7 in terms of x and its derivatives. 
However, it cannot be done straightforwardly because x has singularities at 
values of its argument 7 equal to and 1. Due to momentum conservation 
7(1) = 0, therefore it must be 

X(7(l)) = X(0) = 1 (26) 

which contradicts to its actual behaviour x ~ ols/M (see (|T^). This contra- 
diction can be corrected by regularization: x ~ OLgjM must be replaced at 
small M by 

X~«s/(M + «,). (27) 



Besides the regularization (pTf ), x(^) must be regularized at M ~ 1 . It 
cannot be done with using x(^) at M ~ 1/2 or so because the solution of 
the BFKL is unstable at M = 1/2: the NLO corrections (the second term in 
rhs of (0)) are quite comparable to the leading order contribution. So, in 
such regularization is done through introducing a new parameter A which is a 
"true" value of the Pomeron intercept so that the new improved ^ is supposed 
to be stable at vicinity M = A. Thus, in the context of [0], the expressions 
for new gluon anomalous dimensions contain a presently unknown parameter 
A which must be fixed from experimental data. With such regular izat ions 
of Xi the new anomalous dimensions incorporating NLO contributions to all 
powers in the QCD coupling can be really expressed in terms of x and its 
derivatives. 

The next group of works deals with non-perturbative approaches to the 
Pomeron. In [0 , the contribution of the four-gluon vertices to the Pomeron 



ladder is considered. It is known that such a contribution is negligibly 
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small in the perturbative QCD compared to the three-gluon vertices con- 
tribution. However, ||T^ notes that it is proportional to the trace of the 



energy-momentum tensor in the chiral limit of massless quarks. Therefore, 
the gluon ladder with the four-gluon vertices, apart from a small perturbative 
contribution 

A = /^/"'(M^) (28) 



62 J M6 

to the Pomeron intercept A (in Eq. (pSf ) b = (33 — 2nf)/12TT'^), may con- 
tain a non-perturbative contribution which is proportional to the correlator 
< 9'^{x)9';^{y) > of the energy-momentum tensors. Substituting the spectral 
density p of the correlator in the chiral limit, 

p = (3/327r^)M^ (29) 

into Eq. ( p8|) leads to logarithmic dependency 6 on the upper limit of inte- 
gration Mo in (pHl). Estimating 

= 327r^e,ac/{n} - 1) (30) 

where the value for the vacuum energy e^ac = —{0.2AGeV)'^ is taken from 
the sum rule analysis, they obtain 

A = 0.08. (31) 



The modification of the above approach was made in [jTT|. In essence, it 
comes down to considering Eq. (^) in the same non-perturbative context and 
suggesting to replace in ( PBD by M'^{M'^ + kf)"^ where k'f is an infrared 
cut-off in the transverse momenta space for the exchanged gluons in the 
gluon ladder. Treating this cut-off as a scale for applicability of perturbative 
methods and putting to be equal to the inverse size of the instanton, k"^ = 



{Q.GGeV) , work |TT| concludes that the Pomeron intercept 

A = 0.005. (32) 



Although both previous estimate of |T0| and this estimate are in a good agree- 



ment with experimental data, both works give no prescription what to do 
with the BFKL contribution to the Pomeron intercept, though it is greater 
than the obtained non-perturbative contributions. As NLO corrections to 
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the LO BFKL decrease the value of the LO Pomeron intercept (see [Tl|), 
works |T^, [|ll|] suggest that next corrections may decrease its value down to 



zero and only non-perturbative contributions (^T]),(^2|) would have non-zero 
contributions. However, it is not clear what is the accuracy of the predictions 
(|3TD ,(|32D and what could be corrections to them. 

The works I have discussed contain different improvements of the DGLAP 
and the BFKL to make them more consistent. However, what is really needed 
for the region of small x is a two-dimensional evolution equation combining 
evolutions in x and in so that as were running. Such an equation, the 
infrared evolution equation (IREE), has been obtained recently in work |]12 



for studying the small-x behaviour of the non-singlet contribution /^"^ to the 
structure function Fi. The basic idea for constructing such equations is sim- 
ilar to constructing the renormalization group equations (RGE) but instead 
of evolution with respect to the ultraviolet cut-off, it exploits the infrared 
cut-off evolution: when calculating a structure function, let us introduce the 
infrared cut-off in the transverse momentum space for all integrations over 
virtual particle momenta. With such a regularization, the structure function 
is 11- dependent. Instead of fixing one can evolute the structure function 
with respect to it. The difference between the IREE and the RGE is that 
physical quantities like cross-sections do not depend on the ultraviolet cut-off 
but they have to depend on the infrared one, so within the IREE approach, 
results depend on parameter /x. Work predicts a scaling-like behavior 

(33) 

with the exponent 

a = 0.37. (34) 

This value was obtained when both the leading logarithmic (double-logarithmic) 
and non-leading (single-logarithmic) contributions, including the running 
effects, were taken into consideration. Contrary to the DGLAP prediction 
(^, Eq. ( |5B| ) predicts the power-like small-x behaviour for f'^^. However, 
similarly to the DGLAP, a in Eq. (|33|) depends on n/, Aqcd and on in- 
put fi. The technical difference between the Pomeron and f^^ is that the 
Pomeron gluon ladder is replaced by the quark ladder. Such quark ladder 
is a technically simpler object compared to the Pomeron. Thus, apart from 
physical implications of results obtained in ||T2| , one can regard them as tests 




9 



for checking some of the assumptions made in the works I discussed above. 



Earher, a was calculated in works |13| in the leading, double-logarithmic 



(DL) approximation where was fixed. This DL result for the non-singlet 
structure function at small x is analogous to the leading order BFKL pre- 
diction for the singlet structure function (the Pomeron). Accounting for the 



non-leading contributions in was done in a model-independent way. On 
the other hand, it is analogous to modifications of the BFKL made in works 
Therefore, it is possible and interesting to compare them. Having 
done so, we conclude: 

(i) It is impossible to imitate the running as effects by choosing a reasonable 
scale for fixed as in expressions for intercept of f^^ ■ More precisely, the 
value of the scale for as strongly changes when different kinds of non-leading 
contributions are accounted for. So, concerning results of we think that 
the scale setting would have to be done again if a new portion of non-leading 
contributions to the BFKL were accounted for. 

(ii) Contrary to assumptions made in P, ^, results of [jl2| read that virtual- 



ities of ladder partons can be arguments of as in evolution equations only at 
a; ~ 1. When x is small, the argument of as is more complicated. Account- 
ing for non-leading contributions leads to changing the exponent a in (^) 



but it does not change the fact (obtained in where as was fixed) that 
small x-dependence of f^^ is power-like. Thus, we doubt the DGLAP-like 
a;-dependence obtained in for the Pomeron. 

The small-x behaviour (|33D of f^^ involves a new mass scale /i, with fi^ << 
Q^. Before that, one could think of a Regge-like small-x behaviour f^^ ~ 
[s/Q'^)"' instead of (^). This /i-dependence is due to the fact that the QCD 
perturbative methods cannot be used in the region of too small momenta. 
We think, such a dependence should exist for the Pomeron too. 

I am grateful to M.Greco and S.I.Troyan for useful discussions. 
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